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1. Introduction 

Lattice QCD breaks several global symmetries of the (formal) continuum theory. It follows 
that the corresponding lattice currents are not exactly conserved, and require finite renormalization 
to match to continuum currents. A standard method for determining the renormalization factors 
is to impose Ward-Takahashi identities (WTI) which would hold were the symmetry exact. The 
resulting renormalization factors depend on the WTI chosen, although this dependence vanishes 
with some power of the lattice spacing. 

An important example is the breaking of chiral symmetry when using Wilson fermions (possi- 
bly improved). The current is then the non-singlet axial current, and one determines Za by imposing 
a WTI following from SU {2)i x SU {2)r chiral symmetry (assuming two flavors) 

In this note we address the question of how Za depends upon the choice of WTI, and how 
the dependence enters into results for matrix elements of the renormalized axial current, e.g. fji. 
As noted already, the variation in Za is a discretization error, and thus the question is naturally 
addressed using Wilson chiral perturbation theory (WChPT), i.e. ChPT with discretization errors 
incorporated It turns out, however, that the results on this issue in the literature [jsj, are 
incorrect. Here we present a brief summary of our recent reanalysis |^, explaining what was 
missed in the earlier work, and commenting on the significance of our new results. Our analysis 
controls only discretization errors linear in a, and thus is useful only if either the action or the 
currents (or both) are unimproved. The methodology can, however, be extended to higher order in 
a, and the overall conclusions are applicable more generally. 

We also discuss, in tandem, the corresponding issue for ultra-local lattice vector currents. Here 
there is an exactly conserved lattice current, but it is not ultra-local, and often it is computationally 
simpler to use the ultra-local version. Since the latter current is not conserved, it receives a finite 
renormalization, one that can again be determined by imposing an appropriate WTI. 

2. Currents and Renormalization Conditions 

We consider lattice QCD with Nf = 2 Wilson fermions with lattice spacing a. The ultra-local 
flavor non-singlet vector and axial currents are related to renormalized currents as follows: 

<ren = <Loc ^m,Loc(-^) = V^C-^) 7m V(^) (2-1) 

K,ren = A'^j^^ 4,Loc W = FWrMrS^VW (2-2) 

The Z-factors depend on the action as well as on the WTI used to fix them. For simplicity, we work 
here and below in the massless limit, attained by sending k: ^ jq.. Since we work only at linear 
order in a, we do not encounter the possible Aoki-phase or first-order transition at small masses 
induced by discretization errors of 0{a^) [Q]. This simplification also means that we drop all terms 
proportional to am. ' 

The specific WTI that we use to determine the Z-factors are as follows. For the vector current, 
we impose that the pion have the correct "charge": 

{Ap}KrJ^'{p)) = e''''2E (2.3) 

'in practice, attaining zero quark mass requires an extrapolation, or the use of Sclirodinger-functional boundary 
conditions as an infrared regulator. 
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One must have that p ^Q'm order that £ / (since we are in the massless Umit). The renormal- 
ization factor will, in general, depend on p. 

For the axial current we impose "S^A V" between pion states [|l[ 

|dx£"'-£^^^(7r^(;?)|[A^,,„(jo + ?,f)-Ag,,„(jo-f,f)] A^,,„(j)|7r^(^)) 

= lie^^^%K'{p)%,,,{yW{q)). (2.4) 

Here we need either poxq (or both) to be non-vanishing for the right-hand-side to be non-zero, and 
we choose to take both non-vanishing to avoid infrared divergences. Once we have a renormalized 
vector current (which enters on the right-hand-side) we can use this WTI to determine |Za|. Here 
the result can depend not only on the pion momenta but also on the Euclidean time separating 
the two axial currents (which traces back to the size of the region over which the axial rotation is 
applied). 

Other renormalization conditions are sometimes used in practice, but we have chosen these 
two examples as they can be analyzed using ChPT. We also note that, in practice, one cannot 
consider many choices for the external pion momenta and time separations, due to the associated 
computational cost and the degradation of the signal as \p\ increases. In contrast, ChPT allows one 
to study all momenta in the chiral regime. 



3. Mapping Currents into the Effective Chiral Theory 

The required calculation is now clear: map the lattice currents into the effective chiral theory 



and then evaluate the matrix elements appearing in the WTI ( |2.3| ) and ([2.4[). The mapping is done, 
following in two steps, first from the lattice into the Symanzik effective continuum theory, and 
from there into the chiral effective theory. Since Zy and Za are overall constants to be determined 
at the end, the mapping needs to be done for the bare lattice currents. 

In the first step, using the symmetries of the lattice theory, one finds [p : 

<Loc ^ <Sym,Loc = ^ (^ct + ^^v^vT^tv.ct) +0{a') (3.1) 
1 



^Ji.Loc — ^^,Sym,Loc — ^ ( ^/i,ct + ^^A^/J^rt j +0{a'^) (3.2) 

where the continuum bilinears take their usual forms. This result displays all power-law de- 
pendence on a explicitly (here at linear order); a logarithmic dependence still enters through the 
implicit dependence of cy.A and Zy ^ on g{a). Z-factors are needed in these mappings because the 
currents V^ ^.^ and A^ are conserved while the lattice currents are not. The superscript indicates 
that the Zy ^, while containing perturbative contributions to all orders, do not include 0{a) terms. 
In practice, we know Zy ^ only approximately, but this does not matter since, as will be seen below, 
they cancel once we normalize the currents non-perturbatively. 

The mapping of the currents in ( |3.1[]3.2[ ) into Wilson ChPT is more subtle. One cannot simply 
take the WChPT action (including 0{a) terms) and determine the currents using the Noether pro- 
cedure, as done in Ref. [Q], because the symmetries are broken explicitly. Instead, we have used 
the following two methods. 
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Following continuum ChPT, construct the generating functional, including lattice artefacts 
using standard spurion methods, and obtain the currents by taking derivatives with respect to 



sources. One must treat the 0(1) and 0{a) parts of the currents ( |3.lH3.2| ) separately. 



• Write down the most general currents compatible with the symmetries, and impose appro- 
priate WTI by hand. 

The first approach was used in Ref. [^, and more fully justified in [^. We have checked the 
analysis using a more general method. The second approach is more direct, and we used it to check 
that no subtleties had been overlooked in the first method. We find that both approaches agree. 

The resulting mappings are (in agreement with the results of Ref. |^, and dropping here and 
henceforth a ubiquitous "+0{a^)"): 

<ct+«Cv5v7;v,ct - Vieff = <L0 (l + ^«W^45CSW , (3.3) 

A^^.t+acA^M^i; ^A';;^,ff = A^ LO (^l + ^<3lV45Csw(I+I^)) +AaWKCAd^{T\^-I.'^)), (3.4) 

where ^'^^ lo standard, leading-order ChPT currents, £ contains the pion fields in 

the standard way, d = IWqq, and the Wx are unknown low energy coefficients (LECs) associated 
with the lattice artifacts. The coupling csw is the coefficient of the "clover-term" in the Symanzik 
action, and, like cv.a> is only known approximately. This is now seen to be unimportant, however, 
since it multiplies an unknown LEC. Note that for the vector current the cy term in the Symanzik- 
theory current does not lead to an 0{a) contribution to the current in WChPT, while the ca term in 
the axial current does. 

Putting in the overall factors, the renormalized ultra-local vector and axial currents map into 
WChPT as 

~ ^y* „ A^ ~ — A^ „ (3 5) 

In Ref. [Q], the factors of Zy jZy and Zp^/Z^j^ were set to unity, based on an erroneous argument. In 
fact, these factors can differ from unity at 0{a), as will be seen shortly. 

4. Determining the Renormalization Factors 



We now use the mapped currents of ( |3.5[ ) to evaluate the WTI in WChPT. We work in the 
power-counting in which a ~ 7?^, so that a tree-level calculation suffices to obtain the terms linear 



in a. Evaluating the WTI (2.3) we find the simple result 



ZV=Z0 <,en = <eff (4-1) 

Hence, no additional 0(a) terms are introduced by the current renormahzation. We discuss why 
this is in Ref. [S]; see also Ref. [HI]. 



For the axial current, the WTI (2.4) gives, after some calculation. 



= 1---^(W45Csw + WaCa)za(0, Za(0 = l-cosh[?(|75|-|^|)]exp[-|f||p-^|]. (4.2) 
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In contrast to the vector current we do find a non-vanishing correction of 0{a). This correction 
depends on the separation t between the axial currents and upon the external states, as expected. 
Thus, the renormalized ultra-local axial current maps into WChPT as 



,ren 



4a , X / s 

1 - 72 (W^45CSW + WaCa)Za(0 
J 



V,eff • 



(4.3) 



This is our major new result. 

We conclude that the a dependence of derived using symmetries, is supplemented by 

an additional discretization error resulting from the application of the normalization condition at 
non-zero a. There are three distinct cases (recalling that p,q ^0): 

1. p = q. This is the simplest case to implement practically, and leads to Z\{t) = 0. Thus 
it turns out that, in this case, there are no additional 0{a) terms introduced by the current 
normalization. 

2. p parallel to q. Then, for \t\ ^ 1/1/' — q\, the product of cosh and exponential becomes 1 /2, 
and so za ^ 1 /2. 

3. All other non- vanishing p and q. Here, for \t\^ \/\p — q\, the exponential overwhelms the 
cosh and za — > 1- 

We stress that in both the second and third cases za depends on t for non-asymptotic values of t. 
In any case, our main point is that the implementation of the renormalization condition leads, in 
general, to a non-trivial 0(a) correction to the current. This is what was missed in Ref. 

Another feature of the result ( [4.3[ ) is that it depends on the same unknown coefficients as 



appear in the unnormalized current (3.4), namely the products W45CSW and WaCa. Thus if one were 
to do a fit to multiple physical quantities, incorporating the constraints implied by WChPT so as 
to improve the extrapolation to the continuum limit, the inclusion of the correctly normahzed axial 
current would not increase the number of unknown parameters. 



5. Applications 

Expanding the renormalized axial current we obtain the tree-level result for the pion decay 
constant: 

A,tree = f {\ + (W45CSW + W^A^a) [2 - Za(0]) • (5- 1) 

This gives the form of the discretization errors expected in a lattice calculation of /;;:. We stress 
once again that the result depends on the choice of renormalization condition [through za(0]- 
While expected, it is still striking to see such a dependence explicitly. 

Various comments are in order. First, at this order in ChPT, the continuum result is simply /, 
which is correctly reproduced. Second, a consistency condition on the calculation is that the LECs 
can appear in physical quantities only in certain "physical" combinations and W45CSW + WaCa 
is indeed such a combination. Third, the choice of underlying fermion action enters through the 
values of csw and c\, with the choice of the ultra-local current also affecting c\. If only the action is 
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improved then csw = 0, and one still finds, as expected, an 0(a) term because ca 7^ 0. If both action 
and current are improved, then one finds the expected absence of the 0(a) term — independent of 
the choice of the details of the renormalization condition. Finally, we note that in comparison with 
previous results in the literature, Ref. [^] misses both Wa and za terms, while Ref. [Q] effectively 
assumes that za = 0. 

We have also done a one-loop calculation and find 



jGt,l— loop — f 



a ~ 
1 + -^Wm 



1 



16712/^ 



1 + -^Wa2 



M2 



L45 + ^Wa3 



with W\i being a new unknown LEC while 



(5.2) 



Wax = 4(W45Csw + WaCa)[2-Za(0] "^ai = 4(W45Csw + WaCa)[1 (5-3) 



Note that these two coefficients depend only on the physical combination W45CSW + W^aCa of LEC, 
as expected. 

Our one-loop result correctly reproduces the continuum result of Gasser and Leutwyler 
in the limit a ^ 0. At non-zero lattice spacing, however, there appear additional terms of 0(a), 
0{aM^) and 0{aM^\nM^). Quite generally, the coefficient of the chiral logarithm receives a cor- 
rection in form of the factor [1 + qWai/ f^\, so it not only depends on / and the number of flavors, 
but also on the (non-universal) lattice artifacts encoded in the coefficient Wai An exception 
is the third case discussed in the previous section with Za(0 = 1- Here Wa2 = and the chiral 
logarithm is free of 0(a) corrections. 

Note that the combination L45 of Gasser-Leutwyler coefficients enters the one-loop result in 
form of the lattice spacing dependent combination L^f (a) = L45 + aWA3 //^- In order to obtain the 
physically interesting part L45 one has to extrapolate to the continuum limit. 



We close with a final remark on formula ( p.2[ ). For a consistent result to a given order one 
has to specify a power-counting scheme, and the literature usually distinguishes two regimes: (i) 
the GSM regime with m ~ aA^^^ and (ii) the LCE regime with m ~ a^Agf^^. Here GSM stands 
for generically small quark masses [Q] and LCE for large cut-off effects [^. Equation (5.2) is 
only a partial next-to-leading order (NLO) result for the LCE regime since we ignored the 0{a^) 



corrections in the effective action [10| and the effective currents. The NLO result for the GSM 



regime, however, is obtained from ( |5.2| ) by dropping the corrections proportional to Wa2 and Was, 
which are of next-to-next-to-leading order in this regime. 



6. Conclusions 

We have reconsidered the construction and matching of the vector and axial currents in WChPT. 
The explicit chiral symmetry breaking of Wilson fermions compels us to take two aspects into ac- 
count which are not present in continuum ChPT: 

1. The local lattice currents are not conserved, and in general they do not map onto the con- 
served currents in WChPT. In particular, the WChPT currents are not obtained by the Noether 
procedure, because the currents in the Symanzik theory have 0(a) corrections which are not 
related to the effective action. 
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2. The proper matching of the currents has to take into account the finite renormaUzation of 
the local lattice currents. The same renormaUzation conditions that have been employed 
for the lattice currents must be imposed on the effective currents in WChPT. Depending on 
the particular choice for the renormaUzation conditions the expressions for the renormahzed 
currents differ by terms of 0(a). 

A direct consequence is that WChPT predictions for matrix elements of the currents also depend 
on the renormaUzation condition one has adopted. This is not a flaw but rather reflects the fact that 
the lattice data differs too depending on the particular condition one has chosen. What we find is 
that this dependence enters at 0(a), but that it does not introduce any new LECs. 

Our results emphasize the perhaps rather obvious general point that non-perturbative renor- 
maUzation conditions genetically introduce additional discretization errors. These must be (and in 
practice usually are being) accounted for when extrapolating to the continuum Umit. 
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